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An invariance principle for stochastic heat 
equations with periodic coefficients 

Lu Xu 


Abstract 

In this paper we investigate the asymptotic behaviors of solution u(t , •) 
of a stochastic heat equation with a periodic nonlinear term. Such equa¬ 
tion appears related to the dynamical sine-Gordon model. We consider 
the reversible case and extend the central limit theorem for diffusions in 
finite dimensions presented in in to infinite dimensional settings. Due to 
our results, as t —¥ oo, -j^u (t, •) converges weakly to a centered Gaussian 
variable whose covariance operator is explicitly described. Different from 
the finite dimensional case, the fluctuation in x vanishes in the limit dis¬ 
tribution. Furthermore, we verify the tightness and present an invariance 
principle for {eu(e~ 2 t, -)}te[o,T] as e 4, 0. 

1 Introduction and main results 

A general theory of central limit theorem for additive functionals of re¬ 
versible Markov processes is developed in m, based on a martingale- 
decomposition of the targeted functional. This method is extended to 
non-reversible cases in many references, e.g. nu, m, nu and HU- Com¬ 
bined with Ito’s formula, this method can be used to prove the central 
limit theorem for diffusion processes in finite dimensional space R d with 
both periodic and random coefficients, as illustrated in HI. Chapter 9]. 
In the present paper, we extend this strategy to the case of infinite di¬ 
mensional diffusion processes. 

The model studied in this paper is originally defined in [jj] for the 
purpose of describing the motion of a flexible Brownian string in some 
potential field. The model is expressed by a stochastic PDE on a unit in¬ 
terval [0, 1] with a gradient-type nonlinear term and a Neumann boundary 
condition. Precisely, suppose (D, T, P ) to be a complete probability space 
equipped with a filtration of sub-er-fields {J~t}t>o satisfying the usual con¬ 
ditions in jJJJ Chapter 1], and Wt = W(t, •) to be a cylindrical Brownian 
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motion on L 2 ]0, 1] adapted to Tt. We consider 

dtu(t,x) = i d x u(t,x ) — V x (u(t, x)) + W(t,x), t > 0,x £ (0,1), 

* d x u(t, 0) = d x u(t, 1) = 0, t > 0, 

„ m(0, x) = v(x), x £ [0,1]. 

(1.1) 

In (II.111 . 14'(m) = ^jl4(w) where 14(0 = V(x, •) is a family of C 1 functions 
on R indexed by a; £ [0,1], satisfying 

(1) Vm £ R, V x (u) is Borel-measurable in x\ 

(2) SUPj.gJg^J^gj. {|T4 (m)| + IV4(w)|} < oo; 

(3) Vx £ [0,1], 14 is global Lipschitz continuous with the same Lipschitz 
constant. 

(4) Vx £ [0,1], 14 is periodic in u : V x (u) = 14 (u + 1). 

The motivation of considering such a model lies in several aspects. 
First of all, the study on diffusions in R d already shows the links between 
periodic and random coefficients. Thus our results can provide some hints 
on the more complex case with random coefficients, via the standard tech¬ 
nique of the environment seen from the solution. Secondly, by taking 14 
to be cosine functions m appears closely related to the dynamical sine- 
Gordon equation written as 

d t u = + csin(/3u + 6) + 4 (1.2) 

where c, /3 and 9 are real constants and £ is the standard space-time 
white noise. It is the natural dynamic associated to the usual quantum 
sine-Gordon model, whose solution can be constructed via classical Ito’s 
theory on stochastic PDE when spatial dimension is 1, or via Hairer’s 
theory of regularity structures when spatial dimension is 2, see [8]. From 
a more physical perspective, such dynamic in 2-dimension is first discussed 
in [Tj to study the crystal growth near the roughening transition. It is also 
known that <£3> describes globally neutral gas of interacting charges at 
different temperature /?, see [8j. In the present paper, we restrict our 
discussion to the 1-dimensional case. 

Let S(t) be the semigroup generated by with Neumann boundary 
condition. Then, u(t) = u(t , ■) is called a mild solution to (11.11) if for all 
t > 0, 

u(t) = S(t)v + f S(t — r)[—V.'(u(r,-))]dr + f S(t — r)dW r ■ (1.3) 

Jo Jo 

The regularities and the invariant measure of u(t) are discussed in [6j and 
[7]. Due to [6l Theorem 2.1], u(t) uniquely exists in C[0,1] and forms a 
continuous Markov process. Furthermore, if {w 2 -}x 6 [o,il is a 1-dimensional 
Brownian motion whose initial distribution is the Lebesgue measure on R, 
then the reversible measure of u(t) on C[0,1] is an infinite measure given 
by 

V x (v(x))dx^ fj. w (dv), (1.4) 

where stands for the infinite measure on C[0, 1], induced by w x , see 
[51 Theorem 4.1]. Now we are prepared for presenting our main results. 


/ i{dv) = exp < — 


H 
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Theorem 1.1. Suppose that v is a probability distribution on C[0, 1] being 
absolutely continuous with respect to p in m, and N fT 2 stands for 1- 
dimensional centered Gaussian distribution on R with variance a 2 , where 
a is the constant defined in (f3~7tl below. Then, under initial distribution 
v, for any f £ Cf,(C[0,1]) we have 


lim E v 

E 


1- 

t—> OO 


. \ Vt J 



f 

JR 


f(l-y)N a2 (dy) 


= 0 , 


(1.5) 


where 1 is the constant function on [0,1] defined by l(a:) = 1. 

Remark. Different from the result on diffusions in R d , in Theorem II. II 
the fluctuation in x decays under the time scale t ~ 5 and only the height 
of solution remains in the limit distribution. This kind of decay occurs in 
our model because the Laplacian in the equation controls the growth of the 
fluctuation. More precisely, if one adopts Dirichlet boundary condition in 
EH) instead of Neumann condition, the periodicity turns to be unneces¬ 
sary and the solution u(t ) itself converges weakly to a Brownian bridge 
without any rescaling in time, while it is u(t)/Jt in ESJ. 

By showing the tightness we can prove an invariance principle for u(t). 
Theorem 1.2. Choose T > 0 arbitrarily and fix it. Under an initial dis¬ 
tribution v which is absolutely continuous with respect to p, {eu(e~ 2 t),t £ 
[0, T]} converges weakly to a Gaussian process {oBt • 1 ,i 6 [0, T]} as 
t 0, where Bt is a 1-dimensional Brownian motion on [0, T] and a is 
the constant defined in (1X71) below. 

In Section 2, we study the details of the Kolinogrov generator related 
to (11.11) and present Ito’s formula in Corollary 12.31 In Section 3, we give 
the proofs of Theorem 11.11 and Theorem [Q1 

We end this section by introducing some useful notations. Throughout 
this paper, E is the Banach space C[0,1], and Eq is its subspace consist¬ 
ing of twice continuously differentiable functions ip such that i//(0) = 
ip'f 1) = 0. H is the Hilbert space L 2 [0,1], equipped with the inner prod¬ 
uct (■,•) and the corresponding norm || • ||. Furthermore, let eo = 1 and 
e j( x ) = V2cos(ivjx) for j > 1, then ej £ Eo and {ej}jT 0 forms a complete 
orthonormal system of H. 


2 Reduced process and its Ito’s formula 

Before discussing the generator of EH, we first realize a reduction of 
E. Consider an equivalence relation in E such that ui ~ V 2 if and only if 
ui — V 2 equals to some integer-valued constant function. Define v to be the 
equivalence class of v, and E to be the quotient space E/ ~. From now on 
we identify v with its representative v such that u(0) £ [0,1). A function 
f on E which satisfies that f(v + 1) = f(v) can be automatically regarded 
as a function on E, and vice versa. Furthermore, given a function f on E 
which can be continuously extended to a differentiable function f on H 
with Frechet derivative Df, we simply write Df for Df |^,. 

We explain the benefits of adopting such a reduction. To specify the 
initial condition in EH, we sometimes write u v (t ) instead of u(t). Recall¬ 
ing that 14 is periodic with period 1, the uniqueness of the mild solution 
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shows that u v+1 (t) = u v {t ) + 1 almost surely. This makes it well-defined 
of the reduced process 

u{t) — ( u(t )), Vt € [0, + 00 ), 


where the right-hand side means the equivalence class of u(t) for each 
uj £ fi. This reduction then defines an ij-valued Markov process and 
allows us to convert the infinite measure defined in JH into a finite one. 
Precisely, suppose {iy^.} :Ee [o,i] to be a 1-dimensional Brownian motion 
whose initial distribution is the Lebesgue measure on [0,1), then 

n(dv) = i exp 2 J V x (v(x))dx^ n w (dv), (2.1) 


where n w stands for the probability measure of w' x on E and Z is the 
normalization constant. Apparently, n is reversible for u(t). 

For further discussion on the ergodicity and generator, we consider an 
abstract stochastic differential equation in the Hilbert space H, written 
as 

r dX(t) = [AX{t ) - DV(X(t))]dt + dWt, 

\ X(0) = h £ H, 

where A = ^d x with domain 'D(A) = Eo, V(h) = J V x (h(x))dx for every 
h £ H and DV : H — >■ H is the Frechet derivative. As DV : H —> H is 
bounded and Lipschitz continuous, the mild solution X(t) to (12.211 exists 
uniquely in H, see [3]. 

The probability measure n is the ergodic one. Indeed, denote the 
Hilbert-Schmidt norm by || • ||2 and we know that fg HS'(r) |||dr < 00 for 
all t > 0, thus X{t) is a irreducible strong Feller process, see [16]. Such 
X(t) is regular and can have at most one invariant probability measure, 
see [4] 1.4]. Furthermore, since u(t ) in (11.11) coincides with X(t) when 
h = v € E, u{t) and u{t) are also regular. Consequently, 7r is the unique 
invariant probability measure of u{t). 

Next we discuss the generator of the semigroup of ii{t) on T-L = L 2 ( E , 7r). 
We write the inner product in % by (•, -) n and the norm by || • ||^-. Let Vt 
denote the semigroup by u{t) on T-L. Our aim is to define a Markov pregen¬ 
erator on a relatively simple subset and prove that its closure generates 
Vt- Notice that similar discussion can be found in [2], where they consider 
the dissipativity of a kind of perturbed Kolmogorov operators with finite 
invariant measures. We adopt the ideas in their proofs, however in our 
model 


sup 

t>o 


Tr 


J" 

Jo 


S(r)S* {r)dr 


( 00 

sup < t+y 

00 1 h 


1 — e 


-vr z j z t 


n 2j2 


= OO, 


in contract to [3 Hypothesis l(ii)]. 

Let £a(H) be the linear span of all real and imaginary parts of func¬ 
tions on H of the form h 1 —» e ll ' l ’ h ' 1 where l £ Eo. Suppose £a(E) to be 
the collection of functions in £a(H) such that f(v) = f(v + 1) for all 
v G E, then it can be viewed as a sub-algebra of Cb(E). It is not hard 
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to verify that £a(E) is dense in TL (see e.g. [14]). Consider the following 
Ornstein-Uhlenbeck semigroup on £a(H): 

K t f(h)= [ f(S(t)h + l)M Qt (dl) 

J H 

where Qt = fg S(r)S* (r)dr and Nq t is the centered Gaussian law on H 
whose covariance operator is Qt. Notice that for any fi(h) = we 

have (see ®) 

Htfi=e-^ Qtl ' l) f s * w , (2.3) 

therefore TZ t [£ a(E)\ C £a{E), and IZt = TZt\ £ a (e) can be continuously 
extended to a C° semigroup of contraction operators on TL. Define its 
infinitesimal generator [V (£),£) along Hille-Yosida theorem. Recall that 
/ £ £a(E ) can be continuously extended to a differentiable function on H, 
thus Df is a map from E to H. Ito’s formula suggests that £a(E) C £>(£) 
and for / £ £a(E), 

tf{v) = ( A[Df[v)],v) + \jr [D 2 f(i,)] . (2.4) 

Due to [5] Proposition 1.3.3], £a(E) is a core for L. For / £ £a(E ), let 

tj(y) = tf(v) - (Df(v), DV(v)) (2.5) 

The integration-by-part formula for Wiener measure suggests that 

En\\Df \\ 2 = 2{f, —ICf) n . (2.6) 

From (12.611 . ||A/ — ltf\\ 7 r > > A||/||^, which implies that 1C is 

dissipative, and thus closable. Denote its closure on T-L still by 1C, and the 
domain by V(1C). We prove that 1C generates Vt along the ideas in [2]. 

To this end, we only need to prove the maximal dissipativity, or equiv¬ 
alently, for all A > 0 and g £%, the following equation 

A f-Kf = g 

is solvable in V(1C). The next lemma is similar to [2] Lemma 3.3]. 
Lemma 2.1. Suppose that V £ C 2 (H). For g £ £ A {E) consider the 
equation 

A f-Cf+(Df,DV)=g. (2.7) 

There exists some A > 0 such that the solution to HZ3) exists in Cl (E) fl 
V(C). 

Proof. Let E”[-] = E[-|m( 0) = F] and for A > 0 define function / on E as 

poo 

f(v) 4 / e ~ Xr E ”[s(«(r))]dr. 

Jo 

Let \D 2 V\ oa = sup H \D 2 V(h)\ L ( H ,H) and choose some A > \D 2 V\ 00 . Since 
V £ C 2 (H), we have / £ Cl{E) and 

poo 

(Df{v),h)= / e~ Xr E v [(Dg(u(r)), Y{r))]dr 
Jo 
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where Y(r) £ H is the solution to the following PDE 

dY(r) = AY(r)dr - D 2 V(X(r))(Y(r))dt , Y(0) = h. 

Since 7 Z t f(v) = E” [/ (u(t) + f* S(t - r)[_DF(w(r))]dr^J , we have 

\ (Ktf(v) - /(i>)) = i (lltf(v)-E v [f(u(t))]) + j (e* [/(«(*))]-/(*)) 

converging in 'H to (Df, DV) + A/ — g when t—¥ 0. As C is the generator, 
we know that / G X>(£) and (12.71) holds, so that the proof of Lemma [2.1l is 
completed. Furthermore, as ||Y(r)|| < e r ^ D ^°°||/i||, we have the following 
estimate 

l|£,/ll “ S < 2 ' 8 > 
for every A > \D 2 V\ ao . □ 


Now we are prepared to prove the maximal dissipativity. 
Proposition 2.2. ( T>(K,),K'.) is maximally dissipative on PL. 

Proof. Since we only have V £ Cl(H) in the originally assumption, we 
first choose a sequence V m G C 2 (H) such that \D 2 V rrl \ ao are uniformly 
bounded and 

lim E n \\DVm- DV\\ 2 =0. 

m —>-oo 

Fixing g £ £a(E), from Lemma 12.11 and (12.81) . we are able choose some 
A > 0 large enough such that (12.71) is solvable in Cl (E) D T>(C) with V 
replaced by V m . Denote the solution by / m and pick f m ,n G £a(E ) such 
that both fm t n —> fm and Cf m ,n —> Cfm hold in PL as n —> oo. Such 
sequences exist because £a(E) is a core for C. Combining these with 
(EH) , it suggests that 


lim SUp E-k 11 D fm,ni D fm,ri2 || 

n l > n 2 —too 

—2 lim SUp Eft \( y f'm,ni fm,ri2 ) (E fm,ni Efm,n 25-^^)] 

n l > n 2 —too 

<2||Di/|| 00 liminf j \\fm.ri! - fm,n 2 ||tt [E n ||D/ m , ni - Df m ,n 2 || 2 ] 
n li ri 2 —too I 


thus || Dfm tni — Dfm,n 2 II vanishes in PL. Since the Dirichlet form in m 
is closable in PL (see [II] II.3]), we have linin-^oo E n \\Dfm,n — D/ m || 2 = 0. 
Therefore, 

PL- lim tf m ,n = Cfm - (Df m , DV ), 

n—P oo 

so that fm G V(lt) and LCf m = Cf m - ( Df m , DV). 

Finally, as g - (A/ m - /C/ m ) = (Dfm, DV m - DV), 

lim \\g-(\fm-tfm)\\l< lim {\\Dfm\\oo-E„\\DVm-DV\\ 2 } =0. 

m—Poo m — P oo k 


Recalling that £a(E) is dense in PL, we know that (A — LC)[T>()C)\ is also 
dense. The maximal dissipativity follows from Lumer-Phillips theorem 
(see e.g. Q]5]). □ 
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We have proved that ii(t) generates a strong continuous Markov semi¬ 
group Pt on H whose infinitesimal generator is Furthermore 

£a{E), the linear span of exponential functionals, forms a core of K,. 

Based on Proposition m it is not hard to conclude Ito’s formula for 
u{t). 

Corollary 2.3. For all f £ T>(AC), the following equation holds n-a.s. and 
in T-L. 

u(r))dr+ f (Df(u(r)),dW r )■ (2.9) 

Jo 

Proof. If f € £a(E ), (12.91) follows from the classical ltd’s formula easily. 
For general / £ T>(AC), pick f m £ £a(E ) such that f m —> /, AC/ m —> AC/ 
in H. The corollary then follows from (12.61) and Ito isometry. □ 


/(«(*)) = /(«(0)) + [ ICf( 

Jo 


3 Proofs of Theorems 11.11 and 11.2 

Before illustrating the proof of the central limit theorem, we define two 
Hilbert spaces related to the operator AC. For / £ £a(E) let 

\\f\\l = (-Kf,f)* = ^\\Df\\ 2 . 

Let Hi be the completion of £a{E) under || ■ ||i, which turns to be a 
Hilbert space if all / such that ||/||i = 0 are identified with 0. On the 
other hand, let 

1 = < f £ H; ll/H-i = sup {f,g)n< oo 
l gee A (E), || 9 ||i=i 

Let H-i be the completion of I_i under || ■ ||_i, which also becomes a 
Hilbert space if all / with ||/||-i = 0 are identified with 0. Denote by 
(■,-)i and (-,-)-i the inner products defined by polarization in Hi and 
H-i respectively. For more details about Hi and H- i, see ED Chapter 
2 . 2 ], 

Pick some g> £ Eq and consider a functional V v on E defined as 

V v (v) = i f v(x)ip"(x)dx — f Vx’,(y(x))ip(x)dx. 

2 Jo Jo 

V v is the drift of (11.11) in the direction of tp , meaning that 

(u(t),p) = (u{0),ip) + [ V v {u{r))dr + {W t ,p). (3.1) 

Jo 

Since both p" and are bounded, V v £ H clearly. Furthermore, V<? £ 
£a(E) the integration-by-part formula for Wiener measure shows that 

(V",S), = \ jj,-F,Dg(v))Ti{dv) < i \\p\\E n \\Dg\\ < ^||<p||||s||i. 

Thus £ H-i and ||^||_i < ^\\p\\. 
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Similar as (1231) . consider the resolvent equation for A > 0: 


A ft-tf^ = V v . (3.2) 

Taking inner product with f£ in both sides of (13.211 shows that ||/)f||i < 
(V v , f\)-rr- Since u(t) is reversible under 7r and V v £ H_i we have 

supll/C/m-r — sup \\f\ || i < \\V V || -i < oo. (3.3) 

A>0 A>0 


As the solution lies in D(C), we can consider the Dynkin’s martingale 


Mt(t) = /r(*(*)) - 0)) - /‘/C/r(*(r))dr. 

Jo 

Due to El Chapter 2.6], under condition (13.31) there exists some f v £ Hi 
such that converges to f v in Hi as A 0. Now Corollary 12.31 implies 
that 


limE* 

A 4.0 


f\Dr(u(r)),dW r ) 

Jo 


Um2i||/X’ -/•’Hi = 0. (3.4) 


Furthermore, R^(t) = V ip (u(r))dr — M£(t) vanishes in the following 
sense 

, lim lim K(*) 2 ] = (3-5) 

£—>•00 A —t 


Proof of Theorem M.U Combining (EH) , rm and (1331) yields that as t —> 
00 , the limit distribution of T={u(t), ip) is equal to that of 


^=M*(t) = j\Dr(u(r)) + p, dW r ). 

Now the ergodicity of 7r together with the central limit theorem for con¬ 
tinuous martingales (see e.g. jT8j) shows that under initial distribution 
v < fl, 


lim E v 

t —>00 


E 


( 

V Vt 


To 


I 


f{y)NMdy) 


= 0 


(3.6) 


for all / £ C&(R) and 6 £ R, where aj = E n ^Df v + ip|| 2 . 

Finally, setting <p = ej for j = 1,2,... in (13.61) . we find that f ej (it) = 
—(v, e 3 ), so that a ej = 0 for j > 1. It suggests that CEE holds with 

a 2 = al 0 = limE n \\Dfi° + eof. (3.7) 


In (13.71) . f%° = (A — AC) -1 V e ° is the solution to (13.21) . and its Frechet 
derivative should be understood in the generalized sense, as mentioned 
before Corollary 12.31 □ 

Now fix T > 0 and consider the E-valued process {u^ e \t) = eu(e~ 2 t)} t g[o. 
To prove Theorem Q it is sufficient to verify the tightness. 

Proposition 3.1. Suppose that the law of u( 0) is absolutely continuous 
with respect to p, then as e 4 , 0, the laws ofu' e ' are tight in C([0, T], E). 










Proof. Denote the three terms in the right-hand side of CG2) by X '(t), 
Y(t) and Z(t) respectively. Furthermore, let 

pi 00 

X x (t) 4 X(t) - / X(t, x)dx = V(X(f),e,) • ej . 

Jo 3=1 


Define Y ± (t) and Z x similarly. Due to (11.31) . 

( t ) = e(u(e~ 2 t), e 0 ) + eX ± (e~ 2 t) + eY ± (e~ 2 t) + eZ ± (e~ 2 t). (3.8) 


For the first term in the right-hand side, m yields that 


e(u{e 2 t),e 0 ) 


e{u(0),eo) 


r~ 2t 

— e V e ° (ii(r))dr + e{W e - 2 t ,eo). 

Jo 


Therefore | 111 Theorem 2.32] implies that (e(u(e _2 f),eo)}te[o,T] is tight 
when e 4 , 0. As the heat semigroup is contractive, |eA' x (e _2 t, x)\ < 
esupj.grg ji |m(x)|, which suggests that {eX x (e -2 t)} te [ 0iT j vanishes uni¬ 
formly when e 4 , 0. 

The tightness of the last two terms in (ED) can be derived from the 
following estimates. For all p > 1, there exists a finite constant C > 0 only 
depending on p and { 14 } such that for all ti, t 2 £ [ 0 , 00 ), xi, x 2 £ [ 0 , 1 ], 


E |r x (fi,a;i) - T x (t 2 ,a: 2 )| < C(|ti - f 2 | p + |*i - a; 2 | p ); (3.9) 

E |z x (t 1 ,x 1 ) - X x (t 2 ,a; 2 )| < C(\ti - t 2 \^ + \xi - x 2 \ p ). (3.10) 

We only prove 62) here. For (13.101) . since Z ± are Gaussian variables, 

we only need to prove for p = 1 and the computation is similar. Recall 
, — _ 2 -2 . 

that ej(x) = v2cos(7r jx) so that S(t)ej = e * J ej. First suppose that 
ti = t 2 = t and write 5 X = \x\ — X 2 I, 


12 P 


E\Y- L (t,xi)-Y- L (t,x 2 )\ <Gi 


' 00 r * \ 2.2 1 1 2p 

/ e_7r J sm(njS x )\dr 

j=i J ° 1 1 


Let a = {<4 2 J, the right-hand side is bounded by 


Ci 


a £ 00 1 

E-+ E 4-2 

b'=l J j=a +1 J 


2 p 

aS x 1 

< Ci 

7r 7 r 2 a 


2 P 


< C 2 5 P . 


Next let x\ = x 2 = x in dm write St =t 1 — t 2 and compute similarly, 


\2p 


E\Y- L (ti,x)-Y- L (t 2 ,x)\ <C 3 


E 

j =1 


2 - 2 r 

1 - e _7r J 5t 

n 2j 2 


2 P 


< C 4 4 p . 


Therefore (13.91) holds with some constant depending only on p and 14. 

Due to (13.91) and (13.101) . the last two terms in (13.81) are tight in C([0, T ] x 
[0,1],R). Then they are tight in C([0, T],E) because C([0, T] x [0,1],R) 
and C([0, T],E) have equivalent topologies. □ 
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